Edge excitations of bosonic fractional quantum Hall phases in optical lattices 
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The rapid development of artificial gauge fields in ultracold gases suggests that atomic realization 
of fractional quantum Hall physics will become experimentally practical in the near future. While it 
is known that bosons on lattices can support quantum Hall states, the universal edge excitations that 
provide the most likely experimental probe of the topological order have not been obtained. We find 
that the edge excitations of an interacting boson lattice model are surprisingly sensitive to interedge 
hybridization and edge-bulk mixing for some confining potentals. With properly chosen potentials 
and fluxes, the edge spectrum is surprisingly clear even for small systems with strong lattice effects 
such as bandwidth. Various fractional quantum Hall phases for bosons can be obtained, and the 
phases v — 1/2 and v = 2/3 have the edge spectra predicted by the chiral Luttinger liquid theory. 



Fractional quantum Hall (FQH) phases [l|, Q contain 
a wide variety of interesting physics, including topolog- 
ically degenerate ground states, fractional bulk excita- 
tions, and gapless chiral edge excitations. They arise at 
low temperatures when strong magnetic fields are applied 
to high-quality two-dimensional electron gases with low 
carrier concentration. Ultracold gases of neutral atoms 
are being used to investigate several properties of mate- 
rials which can be hard to control precisely in the solid 
state. As these systems are charge-neutral, an ordinary 
magnetic field cannot be used to create the Lorentz force. 
A synthetic magnetic field can be created by rotation, but 
technical issues appear to limit this approach to lower 
field strengths than are necessary for FQH [3], with the 
exception of a recent experiment with few trapped par- 
ticles Q . Much theoretical work has been done for these 
systems, for a review see including edge spectrum 
calculations @. 

Recently, several theoretical 0-0 and experimen- 
tal @, 3 

proposals have been made for stronger syn- 
thetic magnetic fields for ultracold neutral atoms. All 
of them can be used with optical lattices which enhance 
interaction effects and give a larger energy gap above the 
FQH ground state. Theoretical work on lattice systems 
with an effective magnetic field goes back at least to Hof- 
stader's work [ll| on non-interacting particles; the FQH 
phases are strongly interacting, and S0rensen et al. [7( 
showed that in the low flux limit and strong interactions 
the system could be well described by Laughlin's wave- 
function 12). In subsequent work Hafezi et al. 13j con- 



cluded that this could be extended to larger fluxes per 
unit cell by investigating the topological structure of the 
ground state. 

The goal of this work is to understand practical ex- 
perimental conditions for observation of edge states in 
bosonic lattice FQH systems and compare numerical re- 
sults for edge excitations in hierarchy states to the pre- 
diction of chiral Luttinger liquid theory [14) . Convinc- 
ing observation of bosonic FQH states will depend on 
an experimentally viable probe of the topological order; 



while many multi-particle quantities have been used to 
diagnose the topological state in past theoretical work, 
such as ground state degeneracy, bulk energy gap, wave- 
function overlap, band flatness, band Chern number and 
entanglement spectra, these are not yet experimentally 
accessible. Our focus will be on edge excitations, whose 
"universal" aspects contain information about the topo- 
logical order of the system, although a good understand- 
ing of the "non-universal" effects of the lattice and trap 
is crucial for these excitations to provide a clear signal 
that a FQH phase has been obtained. 

We investigate bosonic FQH phases in the simplest lat- 
tice system, hard-core bosons on a square lattice in a uni- 
form magnetic field with equivalent Landau level filling 
v = N/N$. N is the number of bosons and is the 
number of fluxes in the system, measured in units of the 
magnetic flux quantum $ = hc/q, for particles of charge 
q. The modified Bose-Hubbard Hamiltonian is, 



+ h.c, (1) 



where J is the hopping amplitude, at creates a boson 
on site r = (x,y). We use two different gauges for the 
phases a —(a x ,a y ), Landau gauge a = (a, 0) on cylin- 
ders to keep explicit translational symmetry around their 
circumferences, and symmetric gauge, a = (a/2, a/2) on 
squares to keep explicit Z4 rotational symmetry. The 
flux through a plaquette = a/2ir is defined modulo 1 
and can be expressed as an artificial magnetic field B* = 
ri0$o/ a 2 n, where a is the lattice spacing and h the vector 
normal to the lattice plane. At low flux <C 1, a con- 
tinuum description can be used and the system is effec- 
tively in the flat band limit [IH • We focus on larger fluxes 
where the lattice is important. The magnetic length Ib 
is of the same order as the lattice spacing for the fluxes 
we are interested in, Ib — y/hc/(qB*) = a/ \J2im$ ~ a. 

The spectrum of the Hamiltonian ([1]) is computed with 
exact diagonalization. A section of the system is shown 
schematically in Fig. [Ha). The edge excitations in an 
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FIG. 1. (Color online) (a) Part of the square lattice in the 
symmetric gauge with the phases gained when hopping in the 
direction of the arrows. When hopping in the opposite direc- 
tion, the phases are the complex conjugate of those shown, 
(b) The bulk gap in the v = 1/2 phase on a cylinder as a 
function of flux per plaquette, for N = 4 — 6 particles. The 
large red cross is from a clearly different spectrum, indicating 
that it might not be in a FQH phase. 
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TABLE I. Labeling of the edge spectra for one and two edges: 
2 = {h = l,h = if i ^ 2}, 111 = {h = 3,h =0ifi^ 1} and 
so on. All edge levels with energy E <4hv/L e and positive 
momentum p > are shown. For two edges, the levels with 
p < are equivalent to those with p > 0, with L and R 
exchanged. 



infinite system are gapless, but become gapped in a fi- 
nite system. To clearly see them, it is desirable to have 
a large bulk gap AE^/J = (E^ - E GS )/J, where E% 
is the energy of the lowest bulk excitation and E GS the 
ground state energy in the v phase. The bulk gap on 
cylinders at v — 1/2, in flat infinite wells where edge 
modes do not exist, are shown in Fig. [lib). The ground 
state is non-degenerate and the gaps to all excited states 
are comparable to those reported for a torus Q- To get 
edge excitations in the spectra, more sites outside the 
ground-state droplet need to be added where the edge 
waves can propagate. A trapping potential is then essen- 
tial to confine the condensate. 

The optical trapping potential Vtrap(£j y), provides the 
equivalent of the electric field confinement in semiconduc- 
tors. In this paper, we use harmonic traps Vt ra p{x, y) — 
c x x 2 /a 2 + c y y 2 /a 2 to get one edge on squares and two 
edges on cylinders, with c x — in the periodic direction 
of the cylinder. The velocity of a non-interacting parti- 
cle near the edge is v — VVt ra p(r c )|/(n0/ic/a 2 ), where 
r c is the radius of the ground state droplet. The edge 
excitations of FQH phases form chiral Luttinger liquids. 
In the hydrodynamical approach 14|, the Hamiltonian 
of the edge waves in the Laughlin phases v = 1/m, with 
m = 2, 4, 6, ... for bosons, is 



Hi/ m = 27T- y P-kPkj 

fc>o 



(2) 



k+k' 



wher 



Pk 



with [pk,Pk'] = 
Le 1/2 J d6e ike / h p(9), k = ph/L e with p{6) the one- 
dimensional density along the edge, p E N and L e is 
the length of the edge. This is the U(l) Kac-Moody 
algebra, describing a set of k uncoupled harmonic 
oscillators with energy ^ p l p vk and momentum ^ p l P k, 
where l p G N counts the number of excitations in each 



mode. For a single edge, the degeneracies of the edge 
spectrum are 1,1,2,3,5,7,11,15,22,..., see Tab. Q] for 
a labeling of the different states. On a cylinder, there 
are two edges with momentum in different directions, 
right moving (R) and left moving (L). The degen- 
eracies of this edge spectrum are 1,2,5,10,20,... or 
(1), (1,1), (2, 1,2), (3, 2, 2, 3), (5, 3, 4, 3, 5),... if momen- 
turn is also resolved (Tab. fl}. The wavefunctions in 
the microscopic theory of edge states for the Laughlin 
phases are 



M 2 M 2 



(3) 



i<3 



where P(zi) = J2 P (J2i z iY" ■ This is the form of all zero 
energy wavefunctions without a trap. For a small number 
of particles, the edge excitations consisting of a few single 
particle modes n p = ^2 p l p extends a distance Ar < ti p Ib 
outside the ground state droplet. 

First, we consider the circular harmonic trap on a 
square lattice. With an appropriate construction of the 
Z4 symmetric Hamiltonian, a Fourier transform will turn 
the Hamiltonian into block diagonal form with each block 
corresponding to a certain angular momentum L/h = 
0,1,2,3,... mod 4. These momenta k c — AL/r c , with 
AL = L — L GS , are the same as in the Kac-Moody alge- 
bra k = ph/(2irr c ). The analytic edge spectrum for the 
v = 1/2 phase is shown in Fig. [2fa). In a system with 
finite number of particles N, we only expect excitations 
consisting of n p < N single modes. The energy is propor- 
tional to the angular momentum for the trapped phase, 
indicated by the straight dotted line. With only N = 4 
particles on a 11 x 11 site lattice in a trap c a /J = 0.02, 
we get remarkably good edge spectra. At — 1/5, the 
filling fraction of the ground state is v = 1/2. The de- 
generacies 1, 1, 2, 3, 5, 6, 9, 11, 15, ... in the edge spectrum, 
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agreeing with n p < N for v = 1/2, is clearly visible, 
see Fig. [2fb). A straight line is a good fit to the lowest 
excitation at each momenta. The small energy splitting 
is due to finite size effects, most noticable in the states 
that extends furthest in the trap. Note that all of these 
state appear to be edge states; there are no bulk excita- 
tions for AE 1/2 /J = {E v -Ef, s 2 ) < 0.35, where E v is the 
energy of a state in the v phase, much larger than the 
anticipated gap AEy 2 /J « 0.23. 

Possibly the easiest way to experimentally measure the 
edge spectra is with stimulated two-photon Bragg scat- 
tering [l6j]. The probe light is scattered when the en- 
ergy and momentum difference between the two beams 
is resonant with an excitation level in the FQH system. 
Inserting typical experimental values J = h/rtunnei, with 
a tunneling time T tunn ei =0.2 ms and a = 400 nm gives 
an edge velocity v = AE„/(AL/r c ) ss 0.2 mm/s. Re- 
solving the closest energy excitations at the same angu- 
lar momenta might initially be challenging. However, 
the linear relation v should be well within experimental 
reach with a long pulse duration St ~ 10 ms in the Bragg 
spectroscopy during which the condensate has to remain 
stable. Modulation of the trapping potential is another 
commonly used technique to produce small energy exci- 
tations [17] . Detecting edge excitations could potentially 
also be done by mapping out the in-situ density profile of 
the condensate and an excited state, to observe the prop- 
agating edge mode, or performing the analog to tunneling 
interferometry in condensed matter, either between the 
edges of two condensates trapped next to each other [17] 
or by outcoupling atoms from opposite points of the con- 
densate. 

There are several states [Fig. EJJb)] that do not belong 
to the v = 1/2 edge spectrum. The lowest of those states 
has a v — 2/3 filling factor. The v = 2/3 FQH phase 
has a 3-fold degeneracy on a torus and Moller et al. 18 1 
showed it has a good average overlap with the composite- 
fcrmion wavefunction for < 0.3 on a lattice. Vary- 
ing the magnetic flux slightly, states in the same phase 
change their energy in a similar manner, see Fig. [He) 
for some examples. At small flux per plaquette, our sys- 
tem can become too small for some of the excited states 
of v = 1/2 to exist. Upon decreasing the flux further, 
other phases appear with larger ground state filling fac- 
tors, consistent with some of the phases in the Read- 
Rezayi sequence for bosons v — g/2, where g £ N [19] . 
Whether these phases actually are FQH phases are left 
for future studies. Again, around n$ w 0.3, the edge 
spectrum break down and it is unclear what phases exist 
for ?i0 > 0.3. 

The v = 2/3 hierarchical state [20| , consists of two con- 
densates with comoving edge modes at different radii. 
This spectrum has the same degeneracies as two edges 
on a cylinder 1, (1, 1), (2, 1, 2), (3, 2, 2, 3), but states 
within each set () now have the same angular momen- 
tum, not energy. The edge spectra in Fig. [UJd) is at 



LU 



12 3 
L/h mod 4 



0.35... ... = 

25 = ... 

CSJ 

Lifo.2 ... =..■■--' 

0.15 """ ^ 

0.1 " = 
0.05 ... ^ 

— 

12 3 
L/h mod 4 



-0.05?/ 
-0.1 




A® 

f I 

O ft 'OCOOO o 



v=1/2 
v=2/3 
RR? 



• 


0.4 




(d) 


X 


0.35 


A 




O 


0.3 




0.25 




CO 




llT 0.2 




< 




0.15 




0.1 




0.05 



0.1 



0.2 



0.3 



12 3 
L/h mod 4 



FIG. 2. (Color online) (a) Analytic edge spectrum at v = 1/2 
in a trap as a function of angular momenta. Degenerate lines 
are drawn slightly apart for clarity. The solid lines corre- 
sponds to excitations consisting of four single particle modes 
or less. The dashed dotted lines are the additional excitations 
appearing for additional modes, (b)-(d) Edge excitations for 
N = 4 particles on a 11 x 11 square lattice in a circular har- 
monic trap, v = 1/2 black solid lines, v — 2/3 blue dashed 
lines and the possible Read-Rezayi (RR?) phase, red dashed 
dotted lines, (b) Edge spectrum at = 1/5 as a function of 
angular momentum, (c) Energy gap to the v = 1/2 ground 
state for: ground states (o), 1 states (A), 2 states (x) and 
3 states (•), labeled as in Tab. H as a function of flux per 
plaquette. (d) Edge spectrum at = 1/8 as a function of 
angular momentum. 



= 1/8, where the ground state energy for v = 2/3 
is slightly lower than for v = 1/2. The v — 1/2 edge 
excitations are spaced further apart, but the states with 
degeneracies 1,1,2,3,5 can still be clearly distinguished. 
The two U{\) branches of v = 2/3 have very different ve- 
locitiess, but the structure 1, 2, 5 can easily be seen under 
the dotted line corresponding to the lowest energy for the 
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higher angular momentum states. 

Next, we consider a square lattice on a finite cylin- 
der with a harmonic trap in the non-periodic direction, 
which should be a good approximation to an ellipti- 
cal elongated trap. The cylindrical system in Landau 
gauge has a Zi x symmetry along the circumference of 
the cylinder, when a x = 2tt/L x . With an appropri- 
ate construction of the Hamiltonian, a Fourier trans- 
form will turn the Hamiltonian into block diagonal form 
with each block corresponding to a certain momentum 
k x = p x h/L x — ...,—h/L x ,0,h/L x ,... mod h along the 
circumference of the cylinder. These are the same mo- 
menta as in the Kac-Moody algebra with k x = k. 

The edge spectrum on a cylinder in a 1-D harmonic 
trap c v l J = 0.004, is shown in Fig. _3ja). Three addi- 
tional rows are required outside the ground state droplet 
on both edges to get all edge excitations with \p x \ < N. 
No edge excitations with \p x \ > N are found on any 
cylinder. The edge states 111R1L, 1R111L, 1111R and 
1111L are higher up in the spectrum and not shown. 
The discrepancy from the analytical spectrum, shown in 
Fig. _3b), can be explained by two types of finite size 
effects. The lower-than-expected energy of the high mo- 
mentum single particle modes depends on their overlap 
with the opposite edge, the increase is again from finite 
size effects especially in the furthest extending states. 
Bulk excitations are shown with red lines. The bulk 
gap AE^ 2 /J k 0.12 is again larger than anticipated, 

A£f /2 /J w 0.08 in Fig.^b). 

We briefly mention a few results from other trap po- 
tentials. The low-energy spectra on cylinders does not 
change qualitatively in other potential shapes. How- 
ever, the spectra change completely on squares in cir- 



cular traps Vt ra 



the system shows a FQH edge 



for 1.5 < d < 2.5 with the best spectra at d = 2, with 
AE„ oc c r /J, for reasonable trap strengths. Increasing 
the trapping potential favors denser phases and the tran- 
sitions to phases with higher filling fractions v will hence 
occur at larger fluxes n^. 

Lastly, we discuss which other FQH phases in opti- 
cal lattices that potentially can be detected with our 
method. The other phases in the bosonic Laughlin se- 
quence v — 1/4, 1/6, ... cannot be the ground state with- 
out longer range repulsions in a trap [lijj . In a trap, 
their ground state energies are higher than the ener- 
gies presented here. With decreasing magnetic field, we 
do not observe the higher order hierarchical states at 
v = 3/4,4/5,.... But, we believe this is an artifact of 
the small systems studied. The non-abelian phases in 
the Read-Rezayi sequence v = 1,3/2,2,5/2,... [lj_] are 
believed to be the ground state for bosons in a trap un- 
der certain conditions 



21, 221. The 



3/2 phase is 

the simplest bosonic FQH phase that could be used for 
quantum computing [23j . In the other limit, at large 
flux per plaquette ~ 1/2 (1/3), are the lattice specific 
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FIG. 3. (Color online) (a) Edge spectrum for N = 4 particles 
at v = 1/2 on a 9 x 15 site cylinder in a harmonic trap as a 
function of momenta. Apart from the expected edge spectrum 
(black solid lines), bulk excitations (red dash-dotted lines) 
and edge excitations from translated v — 1/2 phases (blue 
dashed and magenta dashed-dotted lines, where the former 
are aliased to a lower momenta) are mixed in. Two-fold de- 
generate states are marked with longer lines, (b) Analytic 
edge spectrum at v — 1/2 for a large system. Degenerate 
lines are drawn slightly apart for clarity. 



FQH phases believed to appear 15, 1_| 24 1. We consid- 
ered lattices with uniform flux; recent work has shown 
convincingly that non-uniform magnetic fields can cre- 
ate the same FQH phases without Landau levels 25 - 3(| . 
The lessons in this work for trapping potentials and ge- 
ometries apply also to these more complex situations. 
For completeness we mention a few recent proposals for 
detecting quantum Hall physics in ultracold gases other 
than by edge properties [3ll434 

The main result of this work is that hard-core repulsion 
in a small simple lattice system of approximately square 
or circular geometry is sufficient to generate clearly re- 
solved edge excitations for the bosonic FQH states at 
filling v = 1/2 and v = 2/3, provided that the condi- 
tions described above on flux per site and harmonic trap 
strength can be achieved; engineering flat or nearly flat 
bands is unnecessary. In the cylindrical case, the edges 
are strongly interacting with each other and the edge ex- 
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citations are much harder to distinguish, which suggests 
counterintuitively that increasing system size by going 
to an elongated trap may not be necessary or even desir- 
able. Observation of the bosonic FQH states discussed 
here would be a logical first step toward even more excit- 
ing new states that can be studied by similar methods. 
The authors acknowledge helpful conversations with J. 
Guzman, G.-B. Jo, R. S. K. Mong and S. Parameswaran 
and support from the ARO Optical Lattice Emulator pro- 
gram. 
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